Peierls transition as a spatially inhomogeneous gap suppression 
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We propose a model of the Peierls transition (PT) taking into account amplitude fluctuations of 
the charge-density waves and spontaneous thermally activated suppression of the Peierls gap, akin 
to the phase slip process. The activation results in the exponential growth of the normal phase with 
increasing temperature. The model fairly describes the behavior of resistance, thermal expansion, 
Young modulus and specific heat both below and above the PT temperature Tp. The PT appears 
to have a unique nature: it does not comprise Tp as a parameter, and at the same time it has 
features of the 1st order transition. The possible basis for the model is activation of non-interacting 
amplitude solitons perturbing large volumes around them. 

PACS Numbers: 71.45.Lr, 71.30.+h, 68.35.Rh 
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Description of the Peierls transition (PT) in quasi 1- 
dimensional (ID) conductors still remains a controversial 
problem. The widely used mean- field (MF) approach 
works poorly, firstly because of strong ID fluctuations. 
E.g., it predicts the PT temperature Tp much above 
the observed value. Say, for the typical compound TaS3, 
the energy gap 2A is 1600 K, and the MF value of Tp 
should be 2A/3.5 = 460 K, while the experimental value 
is 220 K. The large fluctuations reveal themselves well 
above Tp: the value of the pseudogap is close to the low- 
temperature value [1,2], the threshold non-linear conduc- 
tion [3] indicates the charge-density wave (CDW) state 
within the fluctuating volumes. At present only quali- 
tative attempts to explain these experimental facts are 
undertaken. Though certain success is achieved in fit- 
ting the behavior of different values near Tp [4-6], the 
relations involved are semi-empirical, and their physical 
sense is not quite clear. Another treatment of the PT is 
given by the generalized Ehrenfest relation [7] between 
the specific heat, expansivity, Young modulus anoma- 
lies and stress-induced shift of the transition tempera- 
ture. This relation works nice for some materials, e.g. 
for K0.3M0O3 - the blue bronze [4], but fails for others, 
like for TaS 3 [8]. 

The general approach to the CDWs is to consider them 
as a spatially homogeneous state up to Tp. However, 
a recent theoretical study has demonstrated that ther- 
mal fluctuations of the CDW stress may be very large 
at T « Tp, so that the r.m.s. shift of the chemical po- 
tential level from the middle-gap position appears com- 
parable with A [9], and one can expect temporal local 
suppression of the gap. So, it is natural to suppose that 
the Peierls state in the vicinity of Tp should be consid- 
ered as a mixture of the Peierls phase (the CDW) and 
the state with suppressed gap. Studies of noise have re- 
vealed spontaneous phase-slippage (PS) process [10,11] 
in the vicinity of Tp, which also implies local temporal 
suppression of the Peierl gap. The PS is fairly described 
as a thermally activated process [12-14], so it would be 



intriguing to extend the aproach for the description of 
the PT. 

In the present Letter we propose a model in which the 
fluctuations of A are phenomenologically introduced as 
thermally activated local gap suppression (LGS). With 
increasing T this process gives rise to the activation 
growth of concentration of normal phase as exp(— W/T), 
W ^> A. The approach can be extrapolated above the 
transition temperature. PT appears to be smeared out, 
but at the same time has features of the 1st order transi- 
tion. The model appears to describe fairly the behavior 
of various parameters near Tp. The plausible basics for 
the model is the excitation of amplitude solitons perturb- 
ing considerable volumes around them. 

To introduce the fluctuations, we shall use the follow- 
ing relation giving the frequency of the LGS acts per unit 
volume (see also [14]): 



/ = / a exp(-W/T), 



(1) 



where f a is an attempt frequency. Here the essential 
point is that each LGS act results in a temporal nuclc- 
ation of the normal-state volume vq having certain life- 
time r. Then the fraction of the normal-state volume 
due to the spontaneous LGS process is 



v = v rf. 



(2) 



This fraction grows exponentially with increasing T. At 
high enough temperature v becomes of the order of 1, so 
we should take into account the shrinking of the Peierls- 
state volume. So, instead of (2) we obtain: v = uot/(1 — 
v), or 



1 + VQTf ' 



(3) 



This relation, the principal one for our model, gives 
growth of v from to 1 with increasing T. The growth 
has the form of a step centered at T = W/ ln(v rf a ), at 
which v = 1/2. The step is smeared out by w Tp/W 
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FIG. 1. Fits of a typical R(T) curve for TaS 3 (points) with 
Eq. (4); for the solid line W = 6600 K, v rf a = 3 • 10 13 , 
A = 710 K. The broken line shows a similar fit with 5 = 0. 
The inset shows the corresponding logarithmic derivatives. 

(or l/W in the 1/T-scale). Extrapolation of the LGS 
description above Tp gives a way to treat the entire PT. 
Within this approach the transition consists in the LGS- 
induced destruction of the CDW state. 

Let us compare the resulting relation (3) with the ex- 
periment. We shall refer mostly to TaS3,- a typical CDW 
compound, which is among the widely studied quasi ID 
conductors [15]. The large ratio 2A/Tp for TaS3, as well 
as highly anisotropic structure (the anisotropy of conduc- 
tivity is about 100 at room temperature) argue for the 
strong fluctuations of the CDW order parameter. 

We are beginning with the temperature dependence 
of resistance, the most common curve characterizing the 
PT. Within the present model we should calculate the 
resistance of a mixture of two phases with different re- 
sistivities, p c and p n - We shall consider the Peierls-state 
resistivity p c cx cxp(A/T), and the normal state resis- 
tivity p n = A + BT ', where A and B are constants. To 
calculate the resulting resistivity one should consider a 
complex electric connection of the domains of each phase. 
For simplicity, we take the contribution of each phase to 
the resistivity to be just proportional to the volume frac- 
tion of the phase, as if for connection in series [16]: 

P = vp n + (1 - v)p c . (4) 

Note, that the fluctations are known to contribute to the 
conductivity of the Peierls phase due to thermal depin- 
ning of the CDW [14,17,18]. According to the model 
[14] this contribution is governed by the LGS as well and 
grows as exp(— W/T). Here we shall not distinguish it 
from that of the normal phase. 

Fig. 1 presents an example of a fit of the R{T) curve for 
TaS3 with Eq. (4). The fitting is splendid, but one should 
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FIG. 2. 81/1 vs. T for a TaS3 sample. The background ap- 
proximated with a 2-nd order polynomial is subtracted. The 
solid line gives a fit with Eq. (3) with W = 6500 K. 

note that we have taken B < 0, which is unreasonable 
for a metal. Even if we take B = 0, i.e. p n = const, 
the fitting above Tp becomes considerably worse (see the 
broken line). Below we shall introduce a modification of 
the model above Tp explaining the slower drop of the 
CDW fluctuations. 

Let us now probe the model for the values which com- 
monly characterize thermodynamic transitions, such as 
thermal expansion (TE), Young modulus (Y) and spe- 
cific heat c p (see e.g. [5]). 

To perform TE measurements, we have developed an 
interferometric technique for measurements of needle-like 
samples [19,20]. Fig. 2 gives the temperature dependence 
of the relative length change 51/1 for TaS3 in the vicinity 
of the PT. To exclude the contribution of length hys- 
teresis [19], we present the half-sum of the results ob- 
tained upon heating and cooling the sample. A similar 
curve results if we apply electric field exceeding Ep to 
remove metastability each time before measuring I. Ev- 
idently, the curve presented is close to equilibrium. For 
all our measurements cooling below Tp results in a drop 
of length by about 10~ 5 . This result is quantitatively 
similar to that obtained for K0.3M0O3 in the in-chain 
direction [21]. 

TE at the PT has been discussed in [21]. E.g., it could 
be understood within the anharmonic model, taking into 
account the anharmonic effect of the lattice distortion as- 
sociated with the CDW. Without deepening into details, 
we just assume that the length increase with heating is 
proportional to the fraction of the normal phase. Thus, 
the l(T) step should be described by Eq. (3). The fit 
with W — 6500 K is quite nice (Fig. 2, the solid line). 

As another example we consider the Young modu- 
lus temperature dependence, Y(T). One can expect a 
drop of Y due to the same anharmonic effects. So, with 
T — > Tp from above one can expect a decrease of Y pro- 
portional to the fraction of the CDW volume. However, 



2 




100 150 200 250 



i<r° J — , , , ^ , , , ^ 

4.4 4.6 4.8 5 5.2 5.4 5.6 5.8 

1 3 /T, IC 1 

FIG. 3. Arrhenius plot of the deviation of Y from a linear 
dependence Y(T). The slope of the solid line corresponds to 
W = 5600 K. The inset shows the fit of the whole Y(T) curve. 
The broken line takes into account that for high v (T > Tp) 
1 — v oc cxp( — W/2T). The data are taken from [6]. 

this is not the whole effect. 

Let us recall that depinning of the CDW below Tp 
reduces the Young modulus [4,22,23]. (This effect is 
associated with fast relaxation of the CDW deforma- 
tions which in the pinned state contribute to Y [23].) 
As we mentioned above, with increasing T, the fluctu- 
ations result in the spontaneous depinning of the CDW 
[14,17,18], the concentration of the depinned state grow- 
ing as exp(— W/T) [14], which results in the drop of Y 
with approaching Tp from below. Thus, a dip of Y(T) is 
expected at Tp [24] . The value of the depinning drop of 
Y depends on the particular compound. Being anoma- 
lously strong for TaS3 [22], it is not observed for the blue 
bronze, 5Y/Y < 5-1CT 5 [25,26]. Thus, thcdipinr(T) at 
T — > Tp — should be large for TaS3, and much weaker, 
if any, for the blue bronze. This expectation agrees with 
the experiment: inset to Fig. 3 shows Y(T) for TaS3 from 
[6] . Large drop of Y is seen with T approaching Tp both 
from above and from below, whereas only a small dip of 
Y(T) at T — > Tp from below is observed for the blue 
bronze in the in-chain direction [4,5]. 

Fig. 3 shows the dependence 8Y(T) below Tp in the 
Arrhenius axes (a linear dependence Y(T) is subtracted). 
The slope of the solid line gives the activation energy 
5600 K, giving a good fit nearly up to Tp. To fit the 
whole Y(T) curve we present Y as vY n + v p Y p + v r Y ri 
where the indices n, p and r refer to the normal, pinned 
and relaxed states respectively (v + v p + v r = 1). The 
drop of Y in comparison with the normal state can be 
presented as SY = (Y r — Y n )v r + (Y p — Y n )v p , where 
v r = (1 — v) min(/ r exp(— W/T), 1), and v p = (1 — v — v r ). 
The inset to Fig. 3 shows the fit with W = 6000 K, 
fr = 7.8 • 10 11 , V Tf a = 6 • 10 11 , (Y r - Y n ) = 0.0443, 
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FIG. 4. The temperature derivative of the specific heat 
of TaS 3 . The fit is given by d 2 v/dT 2 (Eq. (3)) with 
W = 10000 K. The data are taken from [8]. 

(Y r — Y p ) = 0.0111 in the normalized units. The fit is 
quite nice, but above Tp the fluctuations fall down slower 
than the fit gives, as it could be expected (recall also 
Fig. 1). 

It is clear from the examples above that the PT con- 
sists in a gradual switching to the state as if having higher 
free energy, and thus looks as a smeared-out 1st order 
transition. The common check for the 1st order transi- 
tion is the latent heat. Because of the smearing out one 
can expect a maximum of specific heat c p . A cusp-like 
feature is clearly seen on the c p (T) curve for the blue 
bronze [5] . Recently a similar feature has been observed 
also for TaS3 [8]. Being very faint, it was detected as a 
zigzag pattern on the derivative dc p /dT. Fig. 4 presents 
the data from [8] together with d 2 v/dT 2 , Eq. (3) (the 
background change of c p is approximated with a straight 
line). It is clear that the form of the feature is at least 
approximately described by our model. Other words, 
the CDW formation is accompanied by a smeared out 
step of latent heat whose width is of the order of Tp/W. 
Note that the values of the latent heat Q w 0.25 R- K 
= 5-10 4 J/m 3 [8] (R is the universal gas constant) and the 
length change 81/1 = 10~ 5 appear to be consistant with 
the Clausius-Clapeyron equation dT p /da = —T p ^j/Q if 
one takes dTp/da ~ 1 K/kbar [15,6] (a is the stress along 
the chains) [27]. 

Above Tp one should bear in mind the small sizes of 
the remnant CDW volumes v c . As soon as they shrink 
down below v n , Eq. (2) is no longer valid, because the 
new normal volume due to the LGS cannot exceed v c . 
With simple assumptions at high enough T one can ob- 
tain (1 - v) oc cxp(-VF/2T); the appropriate fit for Y(T) 
is given with a broken line in the inset to Fig. 3. This con- 
sideration also explains the behavior of the R(T) curve 
above T P (Fig. 1). 

Thus, the LGS model fairly describes the temperature 
evolution of the principle parameters in the vicinity of 
the PT. All the fits proposed have transparent physical 
sense. The main parameter of the model, - the energy 
W, is close to the values for the barrier characterizing 
thermally initiated PS [12-14]. Evidently, LGS is gov- 
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erned by the same process as PS. (Some indications of 
the connection between PS and the PT have been given 
in the early works [10,11,14,28].) So, it would be nat- 
ural to consider excitation of dislocation loops [29] as a 
precursor effect below Tp. Excitation of the dislocation 
loops is being considered as a possible origin of softening 
of solids [30] or similar transitions in liquid helium and 
HTSCs [31]. This approach gives critical expansion and 
proliferation of the loops due to their mutual screening. 
The apparent absence of the critical behavior in our case 
could mean that the CDW excitations practically do not 
interact up to Tp. Note that while for a conventional 
crystal the smallest possible radius of a dislocation loop 
is of the order of the lattice constant, for the CDW it is 
of the order of Such an object (i.e. an amplitude 
soliton) covers a volume ~ £ 3 = where £|| and £j_ 

are the in-chain and the transverse amplitude correlation 
lengths respectively; the soliton can perturb a still higher 
volume, where, say, the conductivity is increased. Thus, 
the condition v = 1/2 can be achieved when the concen- 
tration of the excitations is still <C 1/As, where A is the 
CDW wavelength, and s is the area per chain. 

The concentration of the solitons could be estimated 
as exp(— W/T), where W is the energy of such an 

excitation. Then, v ~ £exp(-W/T). At T = T P we 

have jj exp(— W/T) w 1, and come to the estimate: 



W 



ln(£ 3 /As) 



(5) 



With £ 3 /As = 10 3 we obtain T P = W/7, which can give 
an idea of the low value of Tp in comparison with W. 
A higher ratio W/Tp might be obtained if we take into 
consideration the large wavelengths of the fluctuations of 
the CDW stress along the chains. According to [9] they 
can considerably exceed £|| and, consequently, the LGS 
volumes could appear much larger than £ 3 . 

According to our estimates, the excitations would be- 
gin to turn into dislocation loops, only above Tp. So, 
within the model, the metallic state develops at lower 
temperature than the critical behavior is expected to be- 
gin. As far as we comprehend, the approach proposes a 
new type of phase transition, which does not comprise 
Tp as a parameter. At the same time, the PT in a sense 
resembles a 1st order transition. The model successfully 
works both below and above Tp, though further extrap- 
olation of the approach to higher temperatures requires 
further development of the model. The underlying micro- 
scopic mechanisms of the LGS also needs deeper under- 
standing. Though the model requires further grounds it 
gives a limpid insight into the processes inside the CDW 
near T P . 
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